ON THE EXISTENCE OF NON-CENTRAL WISHART DISTRIBUTIONS 
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Abstract. This paper deals with the existence issue of non-central Wishart distributions 
which is a research topic initiated by Wishart (1928), and with important contributions by 
e.g., Levy (1937), Gindikin (1975), Shanbhag (1988), Peddada & Richards (1991). We present 
a new method involving the theory of affine Markov processes, which reveals joint necessary 
conditions on shape and non-centrality parameter. While Eaton's conjecture concerning the 
necessary range of the shape parameter is confirmed, we also observe that it is not sufficient 
anymore that it only belongs to the Gindikin ensemble, as is in the central case. 



1. Introduction 

The general non-central Wishart distribution T(p,uj;a) on the cone of symmetric pos- 
itive semi-definite d x d matrices is defined (whenever it exists) by its Laplace transform 

C(T(p,u;v))(u) = (det{I + au)y p e- tr{u{I+au ^ lu '\ ueS+, (1.1) 

were p > denotes its shape parameter, a G is the scale parameter and the parameter of 
non-centrality equals oj € 5j\ In the case that oj = 0, F(p; a) := T(p, 0; a) is called the central 
Wishart distribution, which had been introduced in 1928 by Wishart [TTj. In 1937, Levy [5] 
showed that T(p; a) on 5^" is not infinitely divisible for invertible a, which means that for 
some sequence of shape parameters pk I 0, T(pi~;o) cannot exist. Gindikin [5j, Shanbhag [TU] 
and Peddada &; Richards [9j3 subsequently showed that for non-degenerate a, 



(det(J + au)y p (1.2) 

can only be the Laplace transform of a non-trivial probability measure for shape parameters 
p belonging to the Gindikin ensemble 



\ 3 -, j = l,2,...,d-2 



d- 1 
— z — , oo 



Aim of this work is to investigate this fundamental existence issue in the non-central case. 
We shall show: 

Theorem 1.1. Let d G N, p > 0, u, a € . The following hold: 
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(i) Suppose a is invertible. If the right side of (jl.ip is the Laplace transform of a non- 
trivial probability measure^ T(p, co; a) on then p G and rank(w) < 2p + 1. 

(ii) Conversely, suppose any of the following conditions hold: 

(b) p < and rank(w) < 2p. 

Then the right side of is the Laplace transform of a non-trivial probability 

measure T(p,uj;a). 

It should be noted that Theorem is not a full characterization of the existence of non-central 
Wishart distributions, because it leaves open the question, whether distributions T(p, uj;a) 
exist with p G {1/2, . . . , ^-^-} and rank(a;) = 2p + 1. This is only an interesting question for 
d > 3, and rank(w) > 1 as the following two corollaries demonstrate. These are immediate 
conclusions of Theorem 11.11 

Corollary 1.2. Suppose that a is invertible and rank(w) < 1. The following are equivalent: 

(i) The right side of (jl.ip is the Laplace transform of a non-trivial probability measure 
r(p,u; a) on 

(ii) P eA d . 

Another trivial consequence holds in low dimensions. Note that Ai = [0, oo) and A2 = 
[5,00): 

Corollary 1.3. Let d < 2, and suppose a to be invertible. The following are equivalent: 

(i) The right side of (jl.ip is the Laplace transform of a non-trivial probability measure 
T(p, lo; a) on S*j". 

(ii) p G Ad. 

We slightly adapt the notation of the recent article by Letac and Massam [7J, and we are 
recollecting a number of fundamental statements thereof below (see section I2.1|) especially 
what concerns basic properties of non-central Wishart distributions. Concerning their main 
statement as well as that of [H], the following important remark is due: 

Remark 1.4. (i) [TJ Proposition 2.3] claims that A^ fully characterizes the existence of 
non-central Wishart distributions. But this (paradoxically) allows the construction 
of Markovian Feller semigroups on which are non-positive, a mere impossibility. 
That's how we obtain the additional necessary conditions on the rank of u in depen- 
dence of p, which suggests that the characterization of [7J is wrong. On the other 
hand, it is obvious that the existence proof of [7J Proposition 2.3, see also Proposition 
2.1 and the subsequent paragraph] is incomplete, as for p < and rank(o;) > 2p, 
the existence of non-central Wishart distributions T(p,uj;a) is not shown there, 
(ii) [HI Theorem 1] prove the necessity of p G (which had been a conjecture by M.L. 
Eaton) under the premise that rank(w) = 1. However, the method of [Sj, which 
involves the theory of zonal polynomials, relies on the non-negativity of the so-called 
generalized binomial coefficients which may be "difficult to prove" in the case that 
rank(w) > 1, see their concluding remark in Section 4]. In contrast, the present 
paper shows with a much simpler argument that p G A^, for all non-central Wishart 
distributions with nondegenerate scale parameter (see the first part of the proof of 
Theorem 11.11 (i) 



It is easy to see that if a 7^ 0, the triviality of T(p, w; a) is equivalent to p = 0, and lj = (in which case 
we have the point mass at 0). 
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Our method for approaching the existence issue is new. We shall see that p G can be 
proved by utilizing the situation in the central case (restated as Theorem I2.7P and Levy's 
continuity theorem as well as a number of elementary facts such as (1) the behaviour of 
T(p,(jj;a) under the the action of the linear automorphism group of and (2) the charac- 
terizing property of the natural exponential family associated with T(p, u; a), see Proposition 
13.11 Concerning the rank condition, we proceed with an indirect argument: Assuming by con- 
tradiction that p < and rank(o;) > 2p + 1 allows the existence of affine Feller diffusions X 
supported on 5j~ (see [2]) whose transition laws are non-centrally Wishart distributed. The 
contradiction p > is derived by observing that X violates a (geometric relevant) drift 
condition, as established in [2]. The latter is a consequence of the fact that the infinitesimal 
generator of a Markovian Feller semigroup satisfies the strong maximum principle, see [2j 
section 4.4]. 

1.1. Program of the paper. In section [2] we deliver notation and recall known facts about 
the existence of non-central Wishart laws (subsection [24]) and Wishart processes on Sf (sub- 
section [22]). The latter section uses a convenient notation for the Laplace transform, such that 
the distribution of Wishart processes can be easily read off from the characteristic exponents 
of the (affine) process, and which could be easily turned into an existence proof alternative to 
the one of [2]. The presentation of section I2T21 is instructive, and is of relevance for the proof 
of Theorem 11.11 Also, for the sake of completeness, we restate the characterization of the 
central Wishart distributions in terms of the Gindikin ensemble in subsection 12.31 and state 
trivial conclusions when a is degenerate (characterization of existence and infinite divisibil- 
ity). Section presents a proof of Theorem 1 1. 11 In Appendix [A] the relation of our definition 
of non-central Wishart distributions to others in the literature is given. 

2. Notation and preliminary results 

Notation 2.1. Throughout the present article the following notation is relevant: 

• R + is the non-negative real line, and M++ is its interior, 

• Md denotes the set of real d x d matrices, and all symmetric ones therein. 

• I is the unit d x d matrix. 

• St is the cone of symmetric positive semi-definite matrices, and denotes its 
interior, the symmetric positive definite matrices. We denote its boundary Sj" \ S'f 
by dSp 

• ti(A) is the trace of a matrix A G Mj, which introduces a scalar product on via 
(x,y) := tr(xy) for x, y G S d . 

• For k = 1, 2, . . . d— 1, C dS^ determines the (non-convex) cones of d x d matrices 
of rank less or equals k. We note that the dual cone of Dk equals = , hence 
the Laplace transform of a finite measure /i(<i£) supported on is defined by 

C(jm)(u):= [ e-<^V(^), ueS+. 

2.1. Facts on non-central Wishart laws. First we recall the existence and basic properties 
of non-central Wishart distributions: 

Lemma 2.2. Let p € A^, a G S"j" and uj G S"j". We have: 

(i) Suppose w = mm T for m G M. d and set S := a/2. If Y ~ Af(m, X), then X : = 
YY T ~ r(l/2 ,w;a) is supported on D\. 
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(ii) If p < and rank(w) < 2p, then the right side of (jl.ip is the Laplace transform of 
a probability measure supported on Dip- 

(iii) If p > then the right side of (jl.ip is the Laplace transform of a probability 
measure F(p,uj;a) on S^. 

(iv) In particular, if p > and if a is invertible, then the density of T(p,uj;a) exist^ 
and we denote it by F(p,uj,o~,£). 

Proof. Proof of (i) is provided in [7J Proposition 3.2]. (ii) This follows from (i) , by taking 2p 
independent normal random variables on M. d with distribution Af(mi, £) where £ = a/2 and 
oj := m\m[ + • • • + m2 p m^. 

Note that if a G St , our definition of non-central Wishart distribution is related to the 
one of [7] in that F(p, to;a) = / y(p,a^ 1 uja^ 1 ;a), the latter being called "general non-central 
Wishart distribution" in [7] (see Appendix[X]for more detailed information). Hence statement 
is a consequence of (7J p. 1400]. 



iv 



Now for each e > we regularize a and a by setting 

a £ := a + el, a £ := (a + el)~ 1 uj(a + el) . 

Then for each e > 0, we pick X £ , an S~f valued random variable according to [JJ Proposition 
2.3] such that 

X £ ~ Y{p,u);a e )(= -y(p, a £ ;a £ )). 
Letting e — > and using Levy's continuity theorem, we infer that X £ converges in distribution 



to some random variable X ~ T(p,io;a). This settles part (iii) and (iv) 



□ 

2.2. On the Fourier-Laplace transform of Wishart processes. A stochastically con- 
tinuous Markov process (X,¥ x ) x£S + on is called affine, if its Laplace transform is expo- 
nentially affine in the state-variable (see [2]). That is, for all (t,x,u) G R+ x (S^) 2 
E[e- {u ' Xt) \X = x] = e -Ht,u)-m,u),x) ^ ueS + 

holds, where the so-called characteristic exponents <fi and ip satisfy a system of generalized 
Riccati equations, 

4>(t,u) = F(i;(t,u)), 0(O,u)=O, (2.2) 

ip{t, u) = R(tp(t, u)), V(0, u) = u, (2.3) 

and F, R are of a specific Levy-Khintchine form, which is particularly simple in the case of 
Wishart processes (which are pure diffusions; for the original definition in terms of stochastic 
differential equations and particular solutions of these SDEs, see [Jj): 

Definition 2.3. An affine process X = (X,F x ) xeS + is a Wishart process on with pa- 
rameter (p > 0,a G Sj,(3 G Md), if its characteristic exponents (</>, ijj) satisfy the following 
Riccati equations: 

<j)(t,u) = 2p(a,tp(t,u)}, <£(0,u)=0 (2.4) 

ip(t,u) = -2ijj(t,u)aijj(t,u) +ip(t,u)f3 + (3 T ip(t,u), ip(0,u)=u. (2.5) 

Using notation and language from [2j Definition 2.3 and the discussion in section 2.1], the 
parameters satisfy the following 



^For a detailed exposition of the densities, which involves the zonal polynomials, we refer to eq. p. 1400] 
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• a equals the diffusion coefficient of X, 

• b = 2pa equals the constant drift of X, and 

• B(x) = f3x + xf3 T equals the linear drift. Note that the drift enters (|2.5f) as its 
transpose B T (u) = (3 T u + u(3. 

By [21 Theorem 2.4] we have: 

Proposition 2.4. Let a £ St and p > 0. 

(i) If P > ^~T~> then for each (3 E M^, there exists a Wishart process on Sj" with param- 
eters (p, a, f3). 

(ii) Conversely, let X be a Wishart process with parameters (p,a,(3). Then p > 

We further recall the established fact [21 Theorem 2.7, equation (2.22)] that for each x E 
(X, ¥ x ) can be realized as (a weak) solution of the stochastic differential equation 

dX t = yfX t dB t Q + Q T dBjyfX t + (2p Q T Q + /3X t + X t f3 T )dt (2.6) 

subject to Xq = x E Sj", for any Q £ which satisfies Q T Q = a. Here B is a d x d 
standard Brownian motion, and \/X denotes the unique matrix square root on the space of 
positive semi-definite matrices. 

In the following we denote by wf the flow of the vector field fix + xf3 T , that is, 

oj? : S+, wf [x] := f&xeF*. 

Its integral uf : St — > for t > is denoted by 



^ : R+ x 5+ -> 5+, erf (x) = 2 / wf (x)ds. 



Proposition 2.5. Let (X,¥ x ) x&s + be a Wishart process with parameter (p,a,(3). Then the 
characteristic exponents (f>, ijj take the form 

4>(t,u) =plogdet (l + ua^(a)) , (2.7) 



tjj(t,u) = eP Tt (u- 1 +<rf(a)) V*. (2.8) 
Consequently, the Fourier- Laplace transform of X is given by 

i v\ / R \~P - tr( z(l+a?(a)z) 1 u?(x)] 

E x [ e -^ X ')] = (det(/ + af(a)z)J e VV * l W (2.9) 
for all z £ ST + iSd- 



Proof. We first solve the generalized Riccati equations (|2.4p - (|2.5p for initial data u £ St. 
Formula (f!T8|) for ip follows from the fact that jf t a~ l {t) = -a _1 (t)^a(t)o _1 (t), see [U Propo- 
sition III. 4. 2 (ii)]. Formula (|2.7|) follows by some elementary algebraic manipulations using 
the rule f t log(det(a(t)) = tr(a -1 (t)^a(t)), see [H Proposition II.3.3 (i)]. 

Concerning the Fourier-Laplace transform (j2.9|) . we infer directly from their closed- form 
solutions (|2.7p - (|2.8p that 4>{t, u) and ip(t, u) allow for analytic extensions to the complex tube 
+ iSd, which we denote by 4>(t,z) and ip(t,z). Hence using analytic continuation, it can 
be seen that the Fourier-Laplace transform of X is given by 

i v\ / a \~P — tr ( (i +za? (a)] zu>? (x)) 

-< z > x ^ = tdet(I + zo-P(a))) e \S J J, (2.10) 



EBERHARD MAYERHOFER 



In order to obtain (|2.9p from equation (|3.5|) it suffices to observe that for all u, 9 G M d the 
following identity 

u(I + 9u)~ l = (I + u9)~ 1 u (2.11) 

holds, whenever either of both sides is well defined. □ 

Combining Lemma 12.21 and Proposition 12.51 we obtain the following result concerning 
Wishart transition kernels and -densities: 

Proposition 2.6. Let (X,¥ x ) x£S + be a Wishart process with parameters (p,a,f3). Sup- 
pose further that the diffusion parameter a / 0. Then for each (t,x) G ^++ X Sf Xf ~ 

r(p,wf (x); erf (a)). If a G , p > and t > 0, then for all x G Xf has a Lebesgue 
density 

f t (x,0 = F(p,^(x),^(a),0. 

2.3. On the central Wishart case. In the following we restate the characterization of the 
central Wishart laws by using [5]: 

Theorem 2.7. Let d > 2, a G + and p > 0. The following are equivalent: 

(i) Formula (|1.2j) is i/ie Laplace transform of a probability measure T(p,uj;o-) on . 

(ii) p G Arf. 



Proof. (i)=>(ii) This is [9l Theorem 1], as we exclude the point mass at 0. The converse 



direction is a special case of Lemma 12.21 (ii) and (iii) □ 



It is important to note that condition (ii) is not necessary, if a is degenerate. In fact, it is 



easy to prove by use of orthogonal transformations (see |http:/ /arxiv.org/abs/1009.3708 ) 



Corollary 2.8. Let r = rank(<r). The following are equivalent: 

(i) Formula (|1.2p is the Laplace transform of a probability measure T(p,u};a) on . 

(ii) p6 A r . 

As a trivial consequence, one has 

Corollary 2.9. The following are equivalent: 

(i) T(p; a) is infinitely divisible. 

(ii) rank(<r) = 1. 



3. Proof of Theorem 11.11 

Let (p,uj,a) G M ++ x 5j" x such that fi := F(p,uj;a) is a probability measure, that is, 
eq. (jl.ip holds. The domain of its moment generating function is defined as 

D(fi) := {u G S d | £» := / e"^V(^) < oo}, 

which is the maximal domain to which the Laplace transform, originally defined for u G St 
only, can be extended. It is well known that D([i) is a convex (hence connected) set, and 
we also know that S^ C D(fi). Clearly (/ + era) is invertible if and only if the (symmetric) 
matrix (/ + ^/au^/a) is non-degenerate. Using these facts and the defining equation (jl.ip we 
infer that 

£>(/i) := {u G S d | (J + y^u^) G 5++} = -a- 1 + (3.1) 
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and therefore D(p) is even open. Accordingly, the natural exponential family of [i is the 
family of probability measures 3 



f exp«)/x(dg) 

I cm 



v G -a" 1 + S++ 



We start by stating some key properties of Wishart distributions Q: 

Proposition 3.1. (i) Let p >0,w£ S~t. Suppose X is an S~t -valued random variable 
distributed according to F(p,u;I). Let q G S'j" and set cr := q 2 . Then qXq ~ 
T(p,qujq;aj^. In particular, T(p,uj;I) exists if and onlyifT{p,qujq;a) exists. 

(ii) Let p > 0, a G S'j" 1- w G -Sj" suc/i £/ia£ /i := r(p, a;;/) is a probability measure. 
For v = a -1 — I we have that 

exp(ug)/i(dg) 

r(p, cjcjo-; a . 3.2) 

Conversely, if T(p,aoja;a) is a well defined probability measure, so is \i, and (|3.2p 
holds. In particular, we have that the exponential family generated by fj, is a Wishart 
family and equals 

F(p) = {T{p, aua, a) \ a G S++, cr -1 — I € D(ji)}. 

(iii) Let T(p, ujq; o-q) be a probability measure, where o~q G . Then we have 

(a) F(p, tcoo; ctq) is a probability measure for each t > 0. 

(b) If, in addition, ojq is invertible, then T(p,uj\cr) is a probability measure for each 

ueS+,aeS+. 

Proof. Let E be the corresponding expectation operator. By repeated use of the cyclic prop- 
erty of the trace and by the product formula for the determinant, we have 

= det(7 + ou)~ l exp(— tr(uq(I + quq)~ 1 q~ 1 quq)) 
= det(I + au)~ l exp(— tr(u(I + au)~ 1 qujq)), 



which proves assertion (i) Next we show (ii) We note first, that by f)3. 1 1) we have that 
v = a" 1 — 1 G D(fi). Hence exponential tilting is admissible. Furthermore, we have 

( e -^ u+v ^r(p,uj; I)(d£) = det(l + (u + v))~ p exp(- tr((u + v)(l + u + (3.3) 
and setting v = a -1 — lwe obtain 

1 + U + V = <7 _1 (1 + 0~u). 

Hence the first factor on the right side of eq. (|3.3p is proportional to det(l + au)~ p . It remains 
to show that 

- tr((u + v)(l + u + t>)~M = c + tr(«(l + au)~ 1 aua) (3.4) 



^In order to avoid confusions with calculations in the proof of the upcoming proposition, we change here 
from u notation to v, because u denotes the Fourier-Laplace variable in this paper. 

^Some related properties can be found in Letac and Massam [7], but in a different notation. More detailed 
information may be found in Appendix [A"l 

^Expressed in geometric language, we say that the pushfoward of F(p,cj; /) under the map £ i-» q£q equals 
r(p,gwg;cr). 
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for some real constant c, because then the right side of (|3.3p is proportional to the Laplace 
transform of T(p, auja; a). To this end, we do some elementary algebraic manipulations: 

-(tt + v)(I + u + v)~ 1 lj = -(u-l+ cr _1 )[cr _1 (l + au)}" 1 ^ 

= -(-1 + cr" 1 + u)(a~ l + u)~ 1 oj 
= —uj + (<t — a)u + (o - " 1 + u)~ 1 uj 
= (a - I)u - ^((j" 1 + u)^" 1 + uy 1 ^ + (a" 1 + u)~ 
= (a — I)uj — au((j~ l + u)^ 1, 
= (a — I)uj — <ju{I + au)~ ] 



We set now c := tr((cr — I)oS) which is the real number we talked about before. Taking trace 
and performing cyclic permutation inside, we obtain (|3.4|) . and therefore the idendity (|3.2p is 
shown. The assertion concerning the exponential family follows by the very definition of the 
latter. 



We may therefore proceed to (iii) which is proved by repeatedly applying (i) and (ii) Let 
r(p, £Jo;o"o) be a probability measure. Then by (ii) also T(p, a^cjoa^ ;J) is one. Let q\ 
such that q\ = o\ G . We may write T(p, o~q ojqGq 1 , I) = T{p 1 q^ 1 {qiaQ l ojQ(r' {) l qi)q^ l ]I), 
and by applying (i) we obtain the pushforward measure T(p, q\&Q ujqUq q\\ a\). By (ii 
have that T(p,q^ a^ 1 ujqOq 1 q^ 1 ; I) is a probability measure as well, and once again by 



we- 



we infer that for all a G S d + , r(p,aq^ 1 aQ 1 u}QaQ 1 q^ 1 a,a) is a probability. We use this fact 
to prove both parts of the assertion. Without loss of generality we assume that a is non- 
degenerate, because in the case a G dS d we may invoke Levy's continuity theoremQ. Setting 
1/Vil and 



q^a^uJoa^q^ 1 



o"o, we see that 



m a 



holds. For ojq G 51 



a ujo 



which allows to conclude (hi)b 



we choose q\ G St such that 

□ 



Finally, we deliver our proof of Theorem 11.11 

Proof. Let p > such that for some uiq G a G <Sj~ + , the right side of (|1.1|) is the Laplace 
transform of a non-trivial probability measure T(p, ujq; a). By Proposition 13.11 (iii)a we have 
that T(p, u>o/n; a) is a probability measure for each n G N. Letting n — > oo and invoking 
Levy's continuity theorem, we obtain that T(p; a) is a probability measure. But then by the 



characterization of central Wishart laws, Theorem 12.71 (ii) we have that p G A^. 



Let now po G \ [^-,00), and let us assume, by contradiction, that there exist (ujq, a) G 



d ^ u d , rank(o;o) > 2po + 1 such that Y{jpq,ojq\cf) is a probability measure. Pick now 



&d x 

i0\ G Sj such that oj 

d— rank(o;o) 



2_ 



- oj\ + ujq has rank(a/ 
By construction 2p\ = rank(cji), 



) := rank(wi) + rank(ajo) = d, and set p\ : = 
and pi G Ad \ [^-,00). Hence Proposition 



12.21 (ii) implies the existence of a non-central Wishart distribution T(pi, wi, a). Note that 
p* := pq + p\ G Arf \ [^-!-, 00) and that by convolution 

T(p*,u*,a) := r(p ,uo,a)*r(p 1 ,uj 1 ,<j) 



Strictly speaking, Levy's continuity theorem applies to characteristic functions. However, in the Wishart 
case, the right side of (|l.ip can even be extended to even the Fourier-Laplace transform with ease, and by 
preserving its functional form. 
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is a probability measure as well. Since u* is of full rank, we have by Proposition 13.11 (iii)b 
that T(p*,uj;a) is a probability measure for all (to, a) £ (S'j - ) 2 - Hence T(p* ,Lo;ta) exists for 
all (t,u>,a) Gl + x (S+) 2 . 

We proceed by reverse engineering of the results of section 12.21 Pick any a E 5j~ \ {0}. 
For each (t,x) G M + x Sj~ we let pt(x,d^) be the probability measure given by the Laplace 
transform 

/ e- {u ^p t (x,dO = (det(I + 2tau))- pt e t <- uiI+2taurlx ), (3.5) 

(cf. (|2.9p for f3 = 0). By the proof of Proposition 12.51 we know that (p(t, u) := p* log(I + 2tau) 
and tp(t,u) := u(I + 2tau)~ satisfy the system of Riccati equations (|2.4p - (|2,5p with f3 = 
0. From a density argument it follows that the function pt(x,dt;) satisfies the Chapman- 
Kolmogorov equation, hence it is the transition function of a Markov process X on Sj" and 
by construction the Laplace transform is exponentially affine in the state variable x. Hence 
X is an affine process in the sense of [2J, with constant drift parameter b = 2p*a and diffusion 
coefficient a. But b = 2p*a (d — l)a, which contradicts the drift condition formulated in 
Proposition I2.4l|(ii)| Therefore rank(wo) < 2po + 1 and we have proved the first part of the 
theorem. 



The second part of the theorem follows from Lemma 12.21 (ii) and (iii) □ 



Appendix A. Remarks on alternative definitions of the Wishart Distributions 

A number of different notations and definitions of Wishart distributions appear in the 
literature. This paper uses several technical tools which are presented in Letac and Massam's 
work [7j, and therefore we have chosen a notation which is closely related to the latter. 

Letac and Massam use instead of T(p,uj; a) the parameterized family j(p, a; a), where uj is 
replaced by a := a~ 1 (jja~ 1 . Accordingly (jl.ip can be written in the form 

£( 7 (p,a;<7))(u) = (det(I + au)y p e- tv ^ I+au ^ laaa \ ueS+. (A.l) 

Note that this requires a to be invertible. Other authors use densities to define Wishart 
distributions. There are, however, two notable disadvantages of using densities rather than 
the Laplace transform or the characteristic function: 

• A density need not always exist: If a is degenerate, T(p, to; a) is not absolutely 
continuous with respect to the Lebesgue measure on To see this we assume for 
a contradiction that T(p,uj;a) has a Lebesgue density, for some a of rank r < d. 
Let X be an 5j~-valued random variable distributed according to T(p, uj;a). Since 
linear transformations do not affect the property of having a density and since the 
non-central Wishart family is invariant under linear transformations (this is easy to 
check), we may without loss of generality assume that a = diag(0, I r ), where I r is 
the r x r unit matrix. Consider the projection 

7T r : X = (Xjj)i<j : j<d l— ^ Ttr{x) := (xij)l<ij<r- 

A simple algebraic manipulation yields that the Laplace transform of ir r (X) equals 

e -trM^) ; ve s+, 

which is the Laplace transform of the unit mass concentrated at 7r r (uj). But the 
pushforward of a measure with density under a projection must have a density again. 
This yields the desired contradiction. 
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• The density for non-central Wishart distributions is complicated, as it is a series 
expansion in zonal polynomials, see [7]. 
As we do not work with the explicit form of the densities of non-central Wishart distribu- 
tions, we did not need to specify the precise form of the latter in Lemma I2.2l|(iv)[ 



Our use of the Laplace transform throughout the paper is due to the use of affine processes, 
which is a class of Markov processes with the key defining property that their Fourier-Laplace 
transform is of exponentially affine form in the state variable (see, subsection !2.2l in particular 
the defining equation (|2.ip ). One of the important consequences of this property is that the 
affine exponents can be determined by solving a system of ordinary differential equations, the 
so-called Riccati Differential Equations (|2.4p - (|2.5p . rather than the Forward-Kolmogorov- 
PDE for the Laplace transform, 

d t ®{t, u, x) = A$(t, u, x), <3?(0, u, x) = exp(- tr(ux)), 

where A denotes the generator of the affine process. In other words, the last equation can 
be solved by an Ansatz of the form Q(t,u,x) := exp(— 4>(t,u) — tr(ip(t, u)), as it leads to the 
Riccati ODEs and let us avoid solving for complicated parabolic PDEs. For more technical 
details and a complete theory of matrix- variate affine processes we refer the interested reader 
to |2]. 

Even though specifying Wishart distributions by means of their densities is very natural, the 
best known construction is by pushforwards of normal distributions under certain quadratic 
forms: Let ■ ■ ■ De a sequence of Revalued normally distributed random variables 



with mean vectors /ij G M. and covariance matrix E. By using Lemma 12.21 (i) we infer that 
the S'j'-valued random variable 3 := Yli=i £i£7 nas distribution T(p,uj]a), where p = 2k, 

oj = Yh=i Vi^J an d cr = 2S. 

In Gupta and Nagar's notation [6] this reads as follows. We define the d x k matrix 
M = (pi, . . . , pk) an d the d x d matrix O := S~ 1 MM T . Using the matrix-variate normal 
distribution, we have a d x k matrix-valued random variable X := (^i,^2 ; • • • ,£,k) which is 
distributed according to Nd,k(M, £ <8> I), and 

XX T = H ~ Wd(k, X, O), (A.2) 

where k is the shape parameter, £ is the scale parameter, and the parameter of non-centrality 
equals O ([6l Theorem 3.5.1]). As with Letac and Massam's class of generalized non-central 
distributions, this imposes that E must be invertible. Accordingly, the Laplace transform of 
H is given by [6l Theorem 3.5.3] 

E [ e -tr(«s)] = det(J + 2Sn)- fc/2 e- tr(e(/+2&irlSu) . (A.3) 

When E = I then we see that W d (k, E, O = MM T ) equals T(k, MM T ; 2E). Let us now have 
a look at the general case, where E/ / and where O can be any positive semidefinite matrix. 
We observe that the right side of eq. ()A.3p can be also written in the form 

E [ e -tr(«S)] = det (J + 2uS) -fc/2 e -tr(0„S(7+2 U S)-i) ) 

which is the notation found in [9l equation (2)]. To see this, one may use for the first factor 
the multiplicativity of the determinant, while the second factor follows from the identity 

[(I + 2T,u)- 1 T,u}- 1 = (Zu)- 1 ^ + 2Eu) = u' 1 ^ 1 + 21 

as well as 

[uX(I + 2iiE)~ 1 ]~ 1 = S^u" 1 + 21 
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and the symmetry of 0, £ and u. 

We finally show that W d (k, E, 9) = T(k/2, (Qa + ct0)/4; cr). The exponent on the right 
hand side of eq. (|A.3p can be rewritten as 

tr(9(J + au)- 1 au/2) = ^({au)' 1 + cr" 1 )- 1 ©) = tr^w" 1 + a^^Scr /2) 

= tr((u X +a : ) ) 

and since 0, £ are positive semidefinite, also 0cr + a® is positive semidefinite. However, for 
invertible £, the map X i— >■ EX + X£ is injective 0, but not surjective □, in general (unless 
S is a multiple of the unit matrix). Hence the class of non-central Wishart distributions 
used in this paper is strictly larger than the classes in the standard literature, if we insist 
on positive semidefinite non-centrality parameters 0. If we, however, allow (not nessesarily 
positive semidefinite) non-centrality parameters of the form G = 2c~ 1 ix!, where to S S^, 
then all the three mentioned Wishart classes coincide. This definition is also in line with 
the quadratic construction of S, see eq. (|A.2p . Note, however, that [9] imposes a symmetric 
non-centrality parameter, which means that their class of Wishart distributions is strictly 
smaller than ours. Hence Corollary 11.21 comprises a more general situation than [9j Theorem 
2]. 

References 



[2 
[3 
[4 
[5 
[6 
[" 
[8 

[o; 

[10 
[11 



M.-F. Bru, Wishart Processes, J. Theoret. Probab., 4 (1991), pp. 725-751. 

C. Cuchiero, D. Filipovic, E. Mayerhofer, and J. Teichmann, Affine processes on positive semi- 
definite matrices. Ann. Appl. Probab., 21 (2), (2011), pp. 397-463. 

D. Duffie, D. Filipovic, and W. Schachermayer, Affine Processes and Applications in Finance, Ann. 
Appl. Probab., 13 (2003), pp. 984-1053. 

J. Faraut and A. Koranyi, Analysis on symmetric cones, Oxford Mathematical Monographs, The 

Clarendon Press Oxford University Press, New York, 1994. Oxford Science Publications. 

S. G. Gindikin, Invariant generalized functions in homogeneous domains, J. Funct. Anal. Appl. 9 (1975) 

50-52. 

A. K. Gupta and D. K. Nagar, Matrix variate distributions, vol. 104 of Chapman & Hall/CRC Mono- 
graphs and Surveys in Pure and Applied Mathematics, Chapman & Hall/CRC, Boca Raton, FL, 2000. 
G. Letac and H. Massam, The noncentral Wishart as an exponential family, and its moments, J. 
Multivariate Anal., 99 (2008), pp. 1393-1417. 

P. Levy, The arithmetical character of the Wishart distribution, Proc. Cambridge Philos. Soc, 44 (1937). 
S. D. Peddada and D. St. P. Richards, Proof of a conjecture of M. L. Eaton on the characteristic 
function of the Wishart distribution, Ann. Probab., 19 (1991), pp. 868-874. 

D. Shanbhag, The Davidson-Kendall Problem and related results on the structure of the Wishart distri- 
bution, Austr. J. Statist., 30A (1988), pp. 272-280. 

J. Wishart, The generalised product moment distribution in samples from a normal multivariate popu- 
lation, Biometrika, 20A (1928), pp. 32-52. 



Deutsche Bundesbank, Research Centre, Wilhelm-Epstein-Str.14, 60431 Frankfurt am Main, 
Germany 

E-mail address: eberhard.mayerhofer@gmail.com 



This can be seen by diagonalizing E. 
9 lt is well known that any linear transformation on Sj~ is of the form X n> gXg T , where g is a real, not 
necessarily symmetric, d x d matrix. 



